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ABSTRACT 

We present the results of time-integrated searches for astrophysical neutrino 
sources in both the northern and southern skies. Data were collected using 
the partially-completed IccCubc detector in the 40-string configuration recorded 
between 2008 April 5 and 2009 May 20, totaling 375.5 days livetime. An unbinned 
maximum likelihood ratio method is used to search for astrophysical signals. The 
data sample contains 36,900 events: 14,121 from the northern sky, mostly muons 
induced by atmospheric neutrinos and 22,779 from the southern sky, mostly high 
energy atmospheric muons. The analysis includes searches for individual point 
sources and targeted searches for specific stacked source classes and spatially 
extended sources. While this analysis is sensitive to TeV-PeV energy neutrinos 
in the northern sky, it is primarily sensitive to neutrinos with energy greater than 
about 1 PcV in the southern sky. No evidence for a signal is found in any of the 
searches. Limits are set for neutrino fluxes from astrophysical sources over the 
entire sky and compared to predictions. The sensitivity is at least a factor of two 
better than previous searches (depending on declination), with 90% confidence 
level muon neutrino flux upper hmits being between E'^dN/dE ~ 2 — 200 x 
10~^2 TeV cm-2 s~^ in the northern sky and between 3-700 x 10~^^ TeV cm"^ s""^ 
in the southern sky. The stacked source searches provide the best limits to speciflc 
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source classes. The full IceCube detector is expected to improve the sensitivity 
to E""^ sources by another factor of two in the first year of operation. 

Subject headings: astrophysical neutrinos, atmospheric muons and neutrinos 



1. Introduction 

Neutrino astronomy is tightly connected to cosmic ray (CR) and gamma ray astron- 
omy, since neutrinos likely share their origins with these other messengers. With a possi- 
ble exception at the highest observed energies, CRs propagate diffusively losing directional 
information due to magnetic fields, and both CRs and gamma rays at high energies are 
absorbed due to interactions on photon backgrounds. Neutrinos, on the other hand, are 
practically unabsorbed en route and travel directly from cosmological sources to the Earth. 
Neutrinos are therefore fundamental to understanding CR acceleration processes up to the 
highest energies, and the detection of astrophysical neutrino sources could unveil the origins 
of hadronic CR acceleration. Whether or not gamma ray energy spectra above about 10 TeV 
can be accounted for by only Inverse Compton processes is still an ope n question. Some ob- 



servations suggest co ntributions from hadronic acceleration processes (IMorlino et al.l 12009 



Boettcher et al.ll2009l ). Acceleration of CRs is thought to take place in shocks in supernova 
remnants (SNRs) or in jets produced in the vicinity of accretion disks by processes which 
are not fully understood. Black holes in active galactic nuclei (AGN), galactic micro-quasars 
and magnetars, or disruptive phenomena such as collapsing stars or binary mergers leading 
to gamma ray bursts (GRBs), all characterized by relativistic outflows, could also b e pow 



erful accelerators. The canonical model for acceleration of CRs is the Fermi model (IFermi 



19491 ). called first-order Fermi acceleration when applied to non- relativistic shock fronts. 
This model naturally gives a CR energy spectrum with spectral index similar to E~'^ at the 
source. The neutrinos, originating in CR interactions near the sou rce, are expected to fo llow 



a similar energy spectrum. More recently, models such as those in lCaprioli et al.l ( 120101 ) can 
yield significantly harder source spectra. In the framework of these models it is possible to 
account for galactic CR acceleration to energies up to the knee, at about Z x 4 x 10^^ eV, 
where Z is the atomic number of the CR. Extragalactic sources, on the other hand, are 
believed to be responsible for ultra-high energy CRs observed up to about 10^° eV. 

The concept of a neutrino teles cope as a 3-dimensiona l matrix of photomultiplier tubes 



(PMTs) was originally proposed by lMarkov fc Zheleznykhl ( 1l96ll ). These sensors detect the 



Cherenkov light induced by relativistic charged particles passing through a transparent and 
dark medium such as deep water or the Antarctic ice sheet. The depth of these detectors 
helps to filter out the large number of atmospheric muons, making it possible to detect the 
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rarer neutrino events. The direction and energy of particles are reconstructed using the 
arrival time and number of the Cherenkov photons. High energy muon neutrino interactions 
produce muons that can travel many kilometers and are almost coUinear to the neutrinos 
above a few TeV. The first cubic-kilometer neutrino telescope, IceCube, is being completed 
at the South Pole. IceCube has a large target mass. This gives it excellent sensitivity to 
astrophysical neutrinos, enabling it to test many theoretical predictions. 



fl2002f ): Ihearned fc Mannheim! feoooh : 



Gaisser et al 



Re views on neutrino s ource s and telescope s can b e fouri d in lAnchordoaui fc Montarulil 

IChiarusi fc Spuriol ^201^ : iBe c kej loOm : I iLiparil (120061 ) : iBednarek et al.l ( l2005l ): lHalzen fc Hooperl 



Recent results on searches for 



neutri no sources have b een published by IceC ube in the 22-string c onfiguration (lAbbasi et al. 



2009a 



flAmbrosio et al.ll2001f ) 



M), AMANDA-I I flAbbasi et al.ll2009cl ). Super-Kamiokande flThranell2009f ). and MACRO 



This paper is structured as follows: Sec. |2] describes the detector. The data sample and 
cut parameters are discussed in Sec. [3l along with the simulation. In Sec. H] the detector 
performance is characterized for source searches. Sec. describes the unbinned maximum 
likelihood search method, and in Sec. Elthe point-source and stacking searches are discussed. 
After discussing the systematic errors in Sec. [71 the results are presented in Sec. [HI Sec. [9] 
discusses the impact of our results on various possible neutrino emission models, and Sec. [TOl 
offers some conclusions. 



2. Detector and Data Sample 

The IceCube Neutrino Observatory will be composed of a deep array of 86 strings 
holding 5,160 digital optical modules (DOMs) deployed between 1.45 and 2.45 km below 
the surface of the South Pole ice. The strings are typically separated by about 125 m with 
DOMs separated vertically by about 17 m along each string. IceCube construction st arted 



with the first string installed in the 2005-6 austral summer (lAchterberg et al.l l2006al ) and 
will be completed by the end of 2010. Six of the strings in the final detector will use high 
quantum efficiency DOMs and a spacing of about 70 m horizontally and 7 m vertically. Two 
more strings will have standard IceCube DOMs and 7 m vertical spacing but an even smaller 
horizontal spacing of 42 m. These eight strings along with seven neighboring standard strings 
make up DeepCore, designed to enhance the physics performance of IceCube below 1 TeV. 
The observatory also includes a surface array. Ice Top, for extensive air shower measurements 
on the composition and spectrum of CRs. 



Each DOM consists of a 25 cm diameter Hamamatsu PMT (lAbbasi et al 



2010af ). elec- 
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tronics for waveform digitization ( lAbbasi et al.ll2009d[ ). and a spherical, pressure- resistant 
glass housing. A single Cherenkov photon arriving at a DOM and producing a photoelectron 
is defined as a hit. A trigger threshold equivalent to 0.25 of the average single photoelec- 
tron signal is applied to the analog output of the PMT. The waveform of the PMT total 
charge is digitized and sent to the surface if hits are in coincidence with at least one other 
hit in the nearest or next-to-nearest neighboring DOMs within ±1000 ns. Hits that satisfy 
this condition are called local coincidence hits. The waveforms can contain multiple hits. 
The total number of photoelectrons and their arrival times are extracted with an iterative 
Bayesian-based unfolding algorithm. This algorithm uses the template shape representing 
an average hit. 

Forty strings of IceCube were in operation from 2008 April 5 to 2009 May 20. The 
layout of these strings in relation to the final 86-string IceCube configuration is shown in 
Fig. [U Over the entire period the detector ran with an uptime of 92%, yielding 375.5 days of 
total exposure. Deadtime is mainly due to test runs during and after the construction season 
dedicated to calibrating the additional strings and upgrading data acquisition systems. 
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Fig. 1. — Overhead view of the 40-string configuration, along with additional strings that 
will make up the complete IceCube detector. 

IceCube uses a simple multiplicity trigger, requiring local coincidence hits in eight DOMs 
within 5 fis. Once the trigger condition is met, local coincidence hits within a readout win- 
dow ±10 fis are recorded, and overlapping readout windows are merged together. IceCube 
triggers primarily on down-going muons at a rate of about 950 Hz in this (40-string) config- 
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uration. Variations in the trigger rate by about ±10% are due to seasonal changes affecting 
development of CR shower s and muon prod uction in the atmosphere, with higher rates 
during the austral summer (ITilav et al.ll2010l ). 



3. Data and Simulation 
3.1. Data Sample 

Traditional astrophysical neutrino point source searches have used the Ear th to block all 



upward traveling (up-going) particles except muons induced by neutrinos, as in lAbbasi et al. 



( l2009bl ). There remains a background of up-going muons from neutrinos, which are created in 
CR air showers and can penetrate the entire Earth. These atmospheric neutrinos have a softer 
energy spectrum than many expectations for astrophysical neutrinos. The m easurement of 



the atm ospheric neutrino spectrum for the 40-string detector is discussed in lAbbasi et al. 



(l2010bl ). A large number of muons produced in CR showers in the atmosphere and moving 
downward through the detector (down-going) are initially mis-reconstructed as up-going. 
These mask the neutrino events until quality selections are made, leaving only a small residual 
of mis-reconstructed events. 

The down-going region is dominated by atmospheric muons that also have a softer 
spectrum compared to muons induced by astrophysical neutrinos. At present, this large 
background reduces the IceCube sensitivity to neutrino sources in the southern sky in the 
sub-PeV energy region. While veto techniques are in development which will enable larger 
detector configurations to isolate neutrino-induced events starting within the detector, point 
source searches can meanwhile be extended to the down-going region if the softer-spectrum 
atmospheric muon background is reduced by an energy selec tion. This was done for the first 



time using the previous 22-string configuration of IceCube (lAbbasi et al.ll2009af ). extending 
IceCube's field of view to —50° declination. In this paper we extend the field of view to —85° 
declination (the exclusion between —85° and —90° is due to the use of scrambled data for 
background estimation in the analysis, described in Sec. [5]). Down-going muons can also be 
created in showers caused by gamma rays, which point back to their source like neutrinos. 
T he possibility f or Ice Cube to detect PeV gamma ray sources in the southern sky is discussed 



m 



Halzen et al.l ( l2009l ). which concludes that a realistic source could be detected using muons 
in the ice only after 10 years of observing. Gamma ray sources will not be considered further 
here. 

Two processing levels are used to reduce the approximately 3.3 x 10^° triggered events 
down to a suitable sample for analysis (see Tabled]). Random noise at the level of about 
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500 Hz per DOM is mainly due to radioactive decays in the materials in the DOMs. The 
contribution to triggered events by this random noise is highly suppressed by the local 
coincident hit requirement. To further reduce the contribution from noise, only hits within 
a 6 /is time window are used for the reconstructions. This time window is defined as the 
window that contains the most hits during the event. About 5% of down-going muons which 
trigger the detector are initially mis-reconstructed as up-going by the first stages of event 
processing. A persistent background that grows with the size of the detector is CR muons 
(or bundles of muons) from different showers which arrive in coincidence. At trigger level 
they make up about 13% of the events. These coincident muon bundles can mimic the hit 
pattern of good up-going events, confusing a single- muon fit. 

A likelihood-based muon t rack reconstruction is first performed at the South Pole (LI 



filter). The likelihood function (lAhrens et al.ll2004j ) parameterizes the probability of observ- 



ing the geometry and timing of the hits in terms of a muon track's position, zenith angle, 
and azimuth angle. This likelihood is maximized, yielding the best-fit direction and position 
for the muon track. Initial fits are performed using a single photoelectron (SPE) likelihood 
that uses the time of the leading edge of the first photon arriving in each DOM. These re- 
constructions yield robust results used for the first level of background rejection. All events 
that are reconstructed as up-going are kept, while events in the down-going region must 
pass an energy cut that tightens with decreasing zenith angle. Events pass this LI filter at 
an average rate of about 22 Hz and are buffered before transmission via a communications 
satellite using the South Pole Archival and Data Exchange (SPADE) system. 

The processing done in the North includes a broader base of reconstructions compared 
to what is done at the South Pole. Rather than just the simple SPE fit, the multiple 
photoelectron (MPE) fit uses the number of observed photons to describe the expected 
arrival time of the first photon. This first photon is scattered less than an average photon 
when many arrive at the same DOM. The MPE likelihood description uses more available 
information than SPE and improves the tracking resolution as energy increases, and this 
reconstruction is used for the final analysis. The offline processing also provides parameters 
useful for background rejection, reconstructs the muon energy, and estimates the angular 
resolution on an event-by-event basis. Reducing the filtered events to the final sample of 



Table 1. Number of events at each processing level for the 375.5 d of livetime. 



Triggered Events 


3.3 X 10^° 


LI Filtered Events 


8.0 X 10^ 


Events in Final Sample 


36,900 
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this analysis requires cutting on the following parameters: 



Reduced log-Iikelihood: The log-likelihood from the muon track fit divided by the 
number of degrees of freedom, given by the number of DOMs with hits minus five, 
the number of free parameters used to describe the muon. This parameter performed 
poorly on low energy signal events. It was found that low energy efficiency could be 
increased by instead dividing the likelihood by number of DOMs with photon hits 
minus 2.5. Both the standard and modified parameters were used, requiring events to 
pass one selection or the other. This kept the efficiency higher for a broader energy 
range. 

Angular uncertainty, a: An estimate of the uncertainty in the muon track direction. 
The directional likelihood space around the best track solution is sampled and fit to 
a paraboloid. The contour of the p araboloid traces a n error ellipse indicating how 



well the muon direction is localized (INeunhofferl |2006[ ) . The RMS of the major and 



minor axes of the error ellipse is used to define a circular error. This parameter is 
effective both for removing mis-reconstructed events and as an event-by-event angular 
uncertainty estimator. 

Muon energy proxy: The average photon density along the muon track, used as 
a proxy for the muon energy. It is calculated accounting for the distance to DOMs, 
their angular acceptance, and average scattering and absorption properties of photons 
in the ice. The energy loss of a muon moving through the detector scales with the 
muon energy above about 1 TeV when stochastic energy losses due to bremsstrahlung, 
pair production, and photonuclear interactions dominate with respect to ionization 
losses. The energy resolution obtained is of the order of 0.3 in the logio of the muon 
energy (at closest approach to the average hit location) for energies between about 
10 TeV and 100 PeV. The energy of the parent neutrino can be inferred from the 
reconstructed muon energy loss in the detector using simulation. Since the interaction 
vertex is often an unknown distance from the detector, the muon in the detection 
volume has already lost an unknown amount of energy. Fig. [2] shows the distribution 
of this energy parameter versus the true neutrino energy for an E'"^ spectrum. Despite 
the uncertainty on the neutrino energy, for a statistical sample of events this energy 
estimator is a powerful analysis tool because of the wide range over which energies are 
measured. 

Zenith-weighted likelihood ratio: The likelihood ratio between an unbiased muon 
fit and a fit with an event weight according to the known down-going muon zenith dis- 
tribution as a Bayesian prior. Applied to up-going tracks, a high likelihood ratio estab- 
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lishes strong evidence that the event is actually up-going and not a mis-reconstructed 
down-going event. 

• NDir: The number of DOMs with direct photons, defined as arriving within -15 ns 
to +75 ns of the expectation from an unscattered photon emitted from the recon- 
structed muon track at the Cherenkov angle. Scattering of photons in the ice causes 
a loss of directional information and will delay them with respect to the unscattered 
expectation. 

• LDir: The maximum length between direct photons, projected along the best muon 
track solution. 

• Zenith directions of split events: The zenith angles resulting from splitting of an 
event into two parts and reconstructing each part separately. This is done in 2 ways: 
temporally, by using the mean photon arrival time as the split criterion; and geomet- 
rically, by using the plane both perpendicular to the track and containing the average 
hit location as the split criterion. This technique is effective against coincident muon 
bundles mis-reconstructed as single up-going tracks if both sub-events are required to 
be up-going. 




r\ I I I I I I I I I I 1^ 

"^2 3456789 10 
log E^/GeV 

^10 ^ 

Fig. 2. — Distribution of the muon energy proxy (energy loss observed in the detector) versus 
the true neutrino energy for an arbitrary normalization E~'^ flux. 
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In the up-going region, all parameters are used. The zenith-weighted likelihood ratio 
and event splitting are specifically designed to remove down-going atmospheric muon back- 
grounds that have been mis-reconstructed as up-going while the other parameters focus on 
overall track quality. 

In the down-going region, without a veto or earth filter, muons from CR showers over- 
whelm neutrino-induced muons, except possibly at high energies if the neutrino source spec- 
tra are harder than the CR spectrum. The aim of the analysis in this region is therefore to 
select high energy, well-reconstructed events. We use the first three parameters in the list 
above as cut variables, requiring a higher track quality than in the up-going range. Energy 
cuts were introduced in the down-going region to reduce the number of events to a suitable 
size, cutting to achieve a constant number of events per solid angle (which also simplifies 
the background estimation in the analysis). This technique keeps the high energy events 
which are most important for discovery. Cuts were optimized for the best sensitivity using 
a simulated signal of muon neutrinos with spectrum. We checked that the same cuts 
resulted in a nearly optimal sensitivity for a softer spectrum in the up-going region 
where low energy sensitivity is possible and for a harder E~^'^ spectra in the down-going 
region. 

Of the 36,900 events passing all selection criteria, 14,121 are up-going events from the 
northern sky, mostly muons induced by atmospheric neutrinos. Simulations of CR air show- 
ers with CORSIKA show a 2.4±0.8% contamination due to mis-reconstructed down-going 
atmospheric muons. The other 22,779 are down-going events from the southern sky, mostly 
high energy atmospheric muons. An equatorial sky map of these events is given in Fig. |3l 



3.2. Data and Simulation Comparison 



Simulation of neutrinos is used for determining e vent selection and calcula ting upper 
limits. The simulation of neutrinos is based on ANIS ( iGazizov &: Kowalskill2005[ ). Deep in- 
elasti c neutrino-nucleon cross sections use CTEQ5 parton distribution functions ( iLai et al. 
2000l ). Neutrino simulation can be weighted for different fluxes, accounting for the probability 
of each event to occur. In this way, the sar ne simulation sam ple can be u sed to represent at- 
mospheric neutrino models such as Bartol (IBarr et al.ll2004j ) and Honda (IHonda et al.ll2007l ) 
neutrino fluxes from pion and kaon decays (conventional flux). Neutrinos from charmed me 
son decays (prompt flux) have been sim ulated according to a variety of models (IMartin et al. 
20031 : lEnberg et al.ll2008l : lBugaevlll989[ ). Seasonal variations in atmospheric neutrino rates 
are expected to be a maximum of ±4% for neutrinos originating near the polar regions. 
Near the equator, atmospheric variations are much smaller and the variation in the number 
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24h 




Fig. 3. — Equatorial skymap (J2000) of the 36,900 events in the final sample. The galactic 
plane is shown as the solid black curve. The northern sky (positive declinations) is dominated 
by up-going atmospheric neutrino-induced muons, and the southern sky (negative declina- 
tions) is dominated by muons produced in cosmic ray showers in the atmosphere above the 
South Pole. 



of events is expected to be less than ±0.5% (lAckermann fc Bernardinill2005l ). 



Atmospheric muon background is simulated mostly to guide an d verify the event selec- 
tion. Muons from CR air showers were s imulated with C ORSIKA flHeck et al.lll998l ) with 
the SIBYLL hadronic interaction model (lAhn et al.l 120091 ). An October polar atmosphere, 
an average case over the year, is used for the COR SIKA simulation, ignoring the seasonal 
variations of ±10% in event r ates (Tilav et al. 201o[) . Muon propagation through the Earth 
and ice are done using MMC (IChirkin fc Rhodell2004r) . Using measurements of the scattering 
and absorption lengths in ice (lAckermann et al.ll2006r) . a detailed simulation propagates the 
photon signal to each DOM (ILundberg et al.l 120071 ). The simulation of the DOMs includes 
their angular acceptance and electronics. Experimental and simulated data are processed 
and filtered in the same way. 

In Fig. m we show the cosine of zenith and in Fig. [5] the muon energy proxy, reduced 
log-likelihood, and angular uncertainty estimator distributions of all events at trigger level, 
LI filter level, and after final analysis cuts for data and Monte Carlo (MC). In these figures, 
the simulation uses a slightly modified version of the poly-gonato model of the galactic CR 
fiux and composition (jHoerandelll2003l ) . Above the galactic model cutoff at Z x 4 x 10^^ eV, 
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a flux of pure iron is used witli an spectrum. Tliis is done because currently CORSIKA 
cannot propagate elements in the poly-gonato model that are heavier than iron. Moreover, 
the poly-gonato model only accounts for galactic CRs and does not fully account for the 
average na e asured flux above 10^^ eV, even when all nuclei are considered (see Fig. 11 in 
Hoerandell toO^ ). These corrections then reproduce the measured CR spectrum at these 



energies. There is a 23% difference in normalization of data and CR muon events at trigger 
level. This normalization offset largely disappears after quality cuts are made. Generally 
good agreement is achieved at later cut levels. Some remaining disagreements are discussed 
further in Sec. [TJ along with the impact of using different CR primary models. 



4. Detector Performance 

The performance of the detector and the analysis is characterized using the simulation 
described in Sec. 13.21 For an E""^ spectrum of neutrinos, the median angular difference 
between the neutrino and the reconstructed direction of the muon in the northern (southern) 
sky is 0.8° (0.6°). Along with more severe quality selection in the southern sky, the different 
energy distributions in each hemisphere, shown in Fig. |6|, cause the difference in these two 
values. This is because the reconstruction performs better at higher energy due to the 
larger amount of light and longer muon tracks. The cumulative point spread function (PSF) 
is shown in Fig. [7] for two energy ranges and compared with simulation of the complete 
IceCube detector using the same quality selection, as well as the median PSF versus energy 
for the two hemispheres. 

The neutrino effective area A'^ is a useful parameter to determine event rates and the 
performance of a detector for different analyses and fluxes. The expected event rate for a 
given differential flux dN/dE is 

Ne.ents = J dE^Af{E^, ^.) ^^^^^^J ^'^^ , (1) 

and is calculable using simulation. The A'^ represents the size of an equivalent detector if 
it were 100% efficient. Figure [H] shows the A"^ for fluxes of i/^ + z/^ and + u^-, for events at 
final selection level. Neutrinos arriving from the highest declinations must travel through the 
largest column depth and can be absorbed: this accounts for the turnover at high energies 
for nearly vertical up-going muon neutrinos. Tau neutrinos have the advantage that the tau 
secondary can decay back into a tau neutrino before losing much energy. 

Although tau (and electron) neutrino secondaries usually produce nearly-spherical show- 
ers rather than tracks, tau leptons will decay to muons with a 17.7% branching ratio 



- 15 - 




Fig. 4. — Distribution of reconstructed cosine z enith at trigger l evel, LI, and fin al cut level 



for d a ta and simula tion of atmospheric muons (IHoerandell |2003[ ) and neutrinos (IBarr et al. 



2004 : iBugaevI 119891 ) . The true cosine zenith distribution of the muons at trigger level is 
also shown. For the final sample, the deficit of simulation seen near the horizon is discussed 
in Sec. [TJ along with comparisons to other models. The deficit is most likely caused by a 
deficiency in knowledge of the CR composition in the region around the knee. More detailed 
comparisons are discussed in Sec [TJ including comparisons to a variety of models. 



(lAmsler et al.ll2008[ l. At very high energy (above about 1 PeV) a tau will travel far enough 
before decaying that the direction can be reconstructed well, contributing to any extrater- 
restrial signal in the muon channel. For the upper limits quoted in Sec. El we must make 
an assumption on the fiavor ratios at Earth, after oscillations. It is common to assume 
'■ '■ = 1:1:1. This is physically motivated by neutrino production from pion 
decay and the subsequent muon decay, yielding : : = 1:2:0. After stan- 
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Fig. 5. — Distributions of muon energy proxy (top row), reduced log-likelihood (middle row), 
and angular uncertainty estimator (bottom row) for the up-going sample (left column) and 
the down-going sample (right column). Each is shown at trigger level, LI, and final cut level 
for data and simulation of atmospheric muons and neutrinos. In the up-going sample (left 
column), all atmospheric muons are mis- reconstructed, and at final level their remaining 
estimated contribution is about 2.4±0.8%. 



dard oscillations ov er astrophysical baselines, this gives an equal flux of each flavor at Earth 
( lAthar et al.l 120001 ). Under certain astrophysical scenarios, the contribution from muon de- 
cay may be suppressed, leading to an observed flux ratio of '■ '■ = 1 : 1.8 : 1.8 
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Fig. 6. — Energy distribution for an arbitrary normalization E'"^ flux of neutrinos as a 
function of declination for the final event selection. The black contours indicate the 90% 
central containment interval for each declination. 
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Fig. 7. — Cumulative point spread function (angle between neutrino and reconstructed muon 
track) for simulated E'"^ neutrino signal events at the final cut level in the up-going region 
(left). Also shown is the final IceCube configuration. The median of the PSF versus energy 
is shown separately for the northern and southern skies (right). The improvement in the 
southern sky is because of the more restrictive quality cuts. 
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( iKashti fc Waxmaru l2005l ) . or the contribution of tau neutrinos co uld be enhanced by the 



decay of charmed mesons at very high energy (lEnberg et al.ll2009l ). For an E spectrum 



and equal muon and tau neutrino fluxes, the fraction of tau neutrino-induced events is about 
17% for vertically down-going, 10% near the horizon, and 13% for vertically up-going. Be- 
cause the contribution from tau neutrinos is relatively small, assuming only a flux of muon 
neutrinos can be used for convenience and to compare to other published limits. We have 
tabulated hmits on both and the sum + assuming an equal flux of each, while 
in the figures we have specified that we only consider a fiux of muon neutrinos. Limits are 
always reported for the fiux at the surface of the Earth. 
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Fig. 8. — Solid-angle-averaged effective areas at final cut level for astrophysical neutrino 
fiuxes in six declination bands for i/^ + z/^ (left) and Ur + Vt (right), assuming an equal flux 
of neutrinos and anti-neutrinos. 



5. Search Method 



An unbinned maximum likelihood ratio method is used to look for a localized, statisti- 
cally significant excess of events above the background. We also use energy information to 
help separate possible signal from the known backgrounds. 



The method follows that of lBraun et al.l ( 120081 ). The data is modeled as a two component 
mixture of signal and background. A maximum likelihood fit to the data is used to determine 
the relative contribution of each component. Given events in the data set, the probability 
density of the i*^ event is 

(2) 



n 



—)B- 
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where Si and Bi are the signal and background probabihty density functions (PDFs), re- 
spectively. The parameter Ug is the unknown contribution of signal events. 

For an event with reconstructed direction Xi = we model the probability of 

originating from the source at circular 2-dimensional Gaussian, 

1 



A^(^0 = ^e , (3) 



wher e (Tj is the angular uncertainty reconstructed for each event individually ( iNeunhoffer 



20061) and | 

^ i ^ s I 

the space angle difference between the source and reconstructed event. 
The PSF for different ranges of at are in Fig. |9l showing the correlation between the estimated 
angular uncertainty and track reconstruction error. 




'0 0.5 1 1.5 2 2.5 3 



Fig. 9. — Angular deviation between neutrino and reconstructed muon direction A\E' for 
ranges in at, the reconstructed angular uncertainty estimator. Fits of these distributions to 
2-dimensional Gaussians projected into A\E' are also shown. The value of ai is correlated to 
the track reconstruction error. A small fraction of events are not well-represented by the 
Gaussian distribution, but these are the least well-reconstructed events and contribute the 
least to signal detection. 

The energy PDF S{Ei\'y, 6i) describes the probability of obtaining a reconstructed muon 
energy Ei for an event produced by a source of a given neutrino energy spectrum E~'^ at dec- 
lination 6i. We describe the energy distribution using 22 declination bands. Twenty bands, 
spaced evenly by solid angle, cover the down-going range where the energy distributions are 
changing the most due to the energy cuts in the event selection, while two are needed to 
sufficiently describe the up-going events, with the separation at 6 = 15°. We fit the source 
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spectrum with a power law E~'^; 7 is a free parameter. The probabihty of obtaining a re- 
constructed muon energy Ei for an event produced by a source with spectral index 7, for 
spectral indices 1.0 < 7 < 4.0, is determined using simulation. Two examples of these energy 
PDFs are shown in Fig. [TUl 




log^ij(Muon Energy Proxy) log^|j(Muon Energy Proxy) 



Fig. 10. — Probability densities for the muon energy proxy for data as well as simulated 
power-law neutrino spectra. Two declination bands are shown: 0° < S < 15° (left) and 
—40° < 6 < —30° (right), representing two of the declination-dependent energy PDFs used 
in the likelihood analysis. There is an energy cut applied for negative declinations. 

The full signal PDF is given by the product of the spatial and energy PDFs: 

S,=Ar{x,)-£{E,\^,5^). (4) 

The background PDF Bi contains the same terms, describing the angular and energy 
distributions of background events: 

Bi = MAun{x^) ■ £{Ei\Atm, 5i), (5) 

where AfAtm{xi) is the spatial PDF of atmospheric background and £{Ei\Atm,6i) is the 
probability of obtaining Ei from atmospheric backgrounds (neutrinos and muons) at the 
declination of the event. These PDFs are constructed using data and, for the energy term, 
in the same 22 declinations bands as the signal PDF. All non-uniformities in atmospheric 
background event rates caused by the detector acceptance or seasonal variation average out 
in the time-integrated analysis. Therefore MAtm{Si) has a flat expectation in right ascension 
and is only dependent on declination. Because the data are used in this way for background 
estimation, the analysis is restricted from —85° to 85° declination, so that any point source 
signal will still be a small contribution to the total number of events in the same declination 
region. 
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The likelihood of the data is the product of all event probability densities: 

N 



i=l 



The likelihood is then maximized with respect to Ug and 7, giving the best fit values and 
7. The null hypothesis is given by = (7 has no meaning when no signal is present). The 
likelihood ratio test statistic is 



TS = -21o£ 



C{n, = 0) 



(7) 



The significance of the analysis is determined by comparing the TS from the real data 
with the distribution of TS from the null hypothesis (events scrambled in r.a.). We define 
the p-value as the fraction of randomized data sets with higher test statistic values than the 
real data. We evaluate th e median sensitivity and u pper limits at a 90% confidence level 



(CL) using the method of iFeldman fc Cousins! ( Il998l ) and calculate the discovery potential 



as the flux required for 50% of trials with simulated signal to yield a p-value less than 
2.87 X 10^'' (i.e. 5cr significance if expressed as the one-sided tail of a Gaussian distribution). 
The distribution of TS for 10 million trials is shown in Fig. [11] for a fixed point source at 
6 = 45°. Distributions with simulated E'"^ signal events injected are included, as well as 
a distribution with 2 degrees of freedom, which is used to estimate the 5a significance 
threshold for calculating the discovery potential since simulating enough scrambled data sets 
requires a large amount of processing time. 

Although E'"^ sensitivities and limits have become a useful benchmark for comparing 
performance, a wide range of other spectral indices are possible along with cutoffs over a 
wide range of energy. To understand the abihty of the method to detect sources with cut-off 
spectra, typically observed in gamma rays to be in the range 1-10 TeV for galactic sources. 
Fig. [12] shows the discovery potentials for a wide range of exponential cut-offs, demonstrating 
the ability of the method to detect sources with cut-off spectra. Typically, cut-offs observed 
in gamma rays are in the range 1-10 TeV for galactic sources. The likelihood fit is still 
performed using a pure power law. 

The likelihood analysis is not only more sensitive than binned methods, but it can also 
help extract astrophysical information. Fig. [12] shows our ability to reconstruct the spectral 
index for power law neutrino sources at a declination of 6°. The effective area is high for a 
broad range of energies here, and the spectral resolution is best. For each spectrum shown, 
the statistical uncertainty (Icr CL) in the spectral index will be about ±0.3 when enough 
events are present to claim a discovery. Spectral resolution worsens for sources farther from 
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Fig. 11. — Distribution of the test statistic TS for a fixed point source at 5 = 45°. Statistical 
errors are shown for the background-only distribution using 10 million scrambled data sets. 
A distribution with 2 degrees of freedom can be used as an approximation. Also shown 
are the distributions when simulated E'"^ signal events are injected. About 12 events on 
average are needed for a discovery at this declination. 



the horizon, to ±0.4 at both 5 
discovery in each case. 



-45° and 6 = 45° when enough events are present for a 



Stacking multiple sources in neutrino astronomy has been an effective way to enhance 
discovery potential and further constrain astrophysical models (lAchterberg et al.l l2006bl : 
Abbasi et al.ll2009d ). We can consider t he ac cumulated signal from a collection of sources 
using a method similar to lAbbasi et al.l (120061 ). Only a modification to the signal likelihood 
is necessary in order to stack sources, breaking the signal hypothesis into the sum over M 
sources: 



tot 



where is the relative theoretical weight, R-^{'y) is the relative detector acceptance for a 
source with spectral index 7 (assumed to be the same for all stacked sources), and 5/ is the 
signal probability density for the i*'' event, all for the j*^ source. As before, the total signal 
events and collective spectral index 7 are fit parameters. The coefficients depend on 
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Fig. 12. — Flux normalization (left) and number of events at the final level (right) for 
discovery potential versus an exponential cutoff for a differential flux parametrized as 
dN/dE = $° ■ (£'/TeV)~^exp(— .E/i^cutofr)- Curves are shown at three representative de- 
clinations. The likelihood fit is still performed using a pure power law. 
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Fig. 13. — Reconstructed spectral index (la shaded area) versus the number of signal events 
injected for three source spectra: E~^'^, E'"^, and E^^. The sources are pure power-laws at 
a declination of 6°. The stars mark the number of events required for a 5a discovery for each 
spectrum. Systematic errors are not included. 
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our prior theoretical assumptions about the expected neutrino luminosity. They are higher 
for sources that are, on theoretical grounds, expected to be brighter. Tables for R^{j), 
given as the mean number of events from a source with dN/dE oc E~'^ , are calculated using 
simulation. The flexibility built into the method by the relative detector acceptance and 
theoretical weights allows us to use source catalogs covering the whole sky and with large 
variations in source strengths, as well as to directly test model predictions. The improvement 
in discovery potential from stacking sources is illustrated in Fig. [TU showing the discovery 
potential flux versus the number of sources stacked. All sources are at the same declination 
in this example. 
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Fig. 14. — Discovery potential flux versus the number of stacked sources, all with an equal 
E~'^ flux and (for example purposes) at a fixed declination of 45°. Although the total flux 
required increases as more sources are added, owing to a higher effective background, the 
flux per source required can be substantially reduced. 



We would also like to consider sources that are spatially extended (with respect to 
the point spread function). For an example of how important this can be, the significance 
observed by the Milagro experiment in the location of the Fermi source J0634.0+1745 (as- 
sociated with the Geminga pulsar) rises from 3.5cr to 6.3o" by fitting for an extended source 
(lAbdo et al.ll2009l ). The only modification to the method required is to convolve the source 
distribution with the point spread function. Since we model our point spread function as a 
circular 2-dimensional Gaussian distribution, it is easy to also model a source as a circular 
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2-dimensional Gaussian of width o"s. The convolution results in a broader 2-dimensional 
Gaussian of width \/of~+~o^ and the likelihood uses this distribution for the signal spatial 
term. The discovery potential flux for a range of source extensions is shown in Fig. US] 
and compared to the (incorrect) hypothesis of a point source. For a source with true ex- 
tent (Ts = 2°, the point-source hypothesis requires nearly a factor of 2 times more flux for 
discovery compared to the correct extended-source hypothesis. 
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Fig. 15. — Discovery potential flux versus the cr^ of an extended source (distributed as a 
2-dimensional Gaussian) with an E'"^ spectrum at a declination of 25°. The case of a point- 
source hypothesis is compared against the correct extended-source hypothesis matching what 
was used to simulate the signal. 



6. Description of the Five Searches 

We have performed five searches: 

1. a scan for the most significant point source in the entire sky; 

2. a search over an a priori defined list of 39 interesting astrophysical objects; 

3. a stacking search for 16 Milagro TeV gamma ray sources, some seen only in coincidence 
with the Fermi-LAT, and one unconfirmed hot spot (17 total sources); 
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4. a stacking search for 127 local starburst galaxies (IBecker et al.ll2009[ ): 



5. a stacking search for five nearby clusters of galaxies (CGs), testing four different models 
for the CR spatial distribution ( jMurase et al.l 120081 ). 



The analyses and event selection procedure were determined before unblinding the r.a. 
of the data. We require a 5a significance for discovery. Final p- values are calculated for each 
search individually. 



6.1. All-sky Scan 

The first search is a scan for the single most significant point source of neutrinos over 
the declination range —85° to +85°. The maximum likelihood ratio is defined continuously 
over the sky, and we sample it on a grid of 0.1° in r.a. and 0.1° in dec. The size of the grid 
is not important as long as it is small compared to the angular resolution of the detector. 
Using a finer grid increases the computation time with no added benefit. A grid size that is 
comparable to or larger than the angular resolution could miss the location of the peaks in 
the significance map, yielding sub-optimal performance. 



6.2. A Priori Source List 

In order to avoid the large number of effective trials associated with scanning the entire 
sky, we also perform a search for the most significant of 39 a priori selected source candidates, 
given in Table |3l These sources have been selected on the basis of observations in gamma 
rays or astrophysical modeling that predicts neutrino emission. We also added the most 
signifi cant location observ ed in the 22-string IceCube configuration (a post-trial p- value of 



1.3% flAbbasi et al.ll2009bf )) 



6.3. Milagro TeV Source Stacking 



The Milagro Collaboration has reported 16 sources of TeV gamma rays (lAbdo et al. 



20071 ). several only after co rrelating with Ge V gamma rays from the Fermi Gamma-ray 



Space Telescope source list ( lAbdo et al.ll2009l ). These sources are promising candidates for 
detecti on by neutrino telescope s. Particularly interesting are sources in t he complex Cygnus 



region (Beacom &: Kist 



Gonzalez- Garcia et al.l 



ei 



20071 ) and six supernova remnant associations (IHalzen et al.ll2008 



20091), including MGRO J1852+01, a hot spot that falls below the 
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significance threshold of the Milagro Coll. to be claimed as a source. If confirmed as a 
source, MGRO J1852+01 c ould contribute a la rge fraction (about 42%) of the total neutrino 
flux from the SNR sources (IHalzen et al.ll2008l ). For the 40-string configuration of IceCube, 
the model of iHalzen et al.l ( 120081 ) predicts 3.0 neutrino events in 375.5 days, following an 
E~'^'^ spectrum with an exponential cutoff at about 600 TeV. 

We performed a stacking search for 17 sources observed in TeV gamma rays by Milagro 
(adding MGRO J1852+01 to the 16 sources which were found significant by the Milagro 
Coll.) using an equal weight for each source in the likelihood. Assuming that neutrino and 
gamma ray fluxes correlate and using these as weights in the likelihood did not appreciably 
improve the sensitivity in this case. Spatial extensions were used in the search for three 
of the sources where measurements were given (also used in the source simulation for limit 
calculations). The largest source was MGRO J2031+41, reported to have a diameter of 
3.0° ±0.9°. 



6.4. Starburst Galaxy Stacking 



Starburst galaxies have a dense interstellar medium and high star formation rates, par- 



ticularly of high-ma ss stars. This 



eads to both high supernova rates and heating of ambient 



dust. The model of iBecker et al.l (|2009[ ) associates the far infrared (FIR) emission with this 



hot dust and the radio emission with synchrotron losses of CR electrons, presumably ac- 
celerated along with hadronic CRs in the elevated number of SNRs. The observed strong 
correlation between the FIR and radio emission points to the high star formation rate as 
the single underlying cause, and should also correlate with the neutrino flux. The increased 
production of CRs and high density of target material are ideal conditions for neutrino pro- 
duction. The starb urst galaxies M82 and NGC 253 have been observed in gam ma rays at 
GeV-TeV energies (lAbdo et al.ll2010l : lAcciari et al.ll2009l : lAharonian et al.ll2005l ) and are the 
only observed steady extragalactic TeV gamma ray sources not associated with AGNs. 

We performed a stacking search for 127 starburst galaxi es, weighting the so urces by 
their observed FIR flux at 60 /im, as compiled in Table A.l in iBecker et al.l (120091 ). 



6.5. Galaxy Cluster Stacking 

Clusters of galaxies (CGs) are another potential source of high energy protons and, 
through interactions with intracluster material (ICM), neutrinos. CGs are the largest gravi- 
tationally bound objects in the Universe and continue to grow through merging and accretion 
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of dark matter and baryonic gas, generating shock fronts on megaparsec scales. The possi- 
bihty for CGs to be s ources of ult r a high energy CRs above 3 x 10^^ eV is described in, e.g.. 



Norman et al.l (Il995[ ): iKang et al.l ( 119971 ). iMurase et al.l (120081 ) discuss the possibility of CGs 
being a significant contribution to the CR spectrum between the second knee at about 3x 10^'' 
eV and the ankle at about 3 x 10^^ eV. They give predictions for neutrinos from five nearby 
{z < 0.03) CGs: Virgo, Centaurus, Perseus, Coma, an d Ophiuchus. Information on l ocation, 
distance, and size of CGs (virial radii) was taken from iReiprich fc Boehringerl ( 12002| ). These 
nearby CGs appear to us as spatially extended objects with virial radii subtending 1.3°-6.9°, 
so an extended spatial distribution of neutrino s is possible. Whereas the distribution of the 
ICM is well known from X-ray observations (jPfrommer &: EnsslinI 120041 ). the distribution 
of CRs is highly uncertain. The distribution of neutrinos is given by the product of the 
CR and IC M distributions . Four CR models have been conside red for neutrino production 
discussed in iMurase et al.l ( 120081 ) and references therein (e.g., ( IColafrancesco fc Blasilll998 



Berezinsky et al. 19971 )): 



Model A: CRs are uniformly distributed within the cluster shock radius, taken to be 
0.56 of the virial radius for the dynamical parameters considered. 

Model B: CRs are uniformly distributed within the virial radius, yielding the most 
conservative neutrino flux distributed over the largest area. 

Isobaric: CRs follow the distribution of thermal gas. 

Central AGN: In a two-step acceleration scenario CRs are accelerated in the central 
AGN up to a maximum energy before diffusing throughout the cluster and possibly 
undergoing further acceleration. For the purposes of IceCube sea rches, this model ca n 



be treated as a point source. This model is discussed in detail by iKotera et al.l ( l2009l ). 



Signal neutrinos were simulated according to each of the four models. We modeled the 
source extensions in the likelihood as 2-dimensional Gaussian distributions with the width 
for each source and each model determined by optimizing for the best discovery potential. 
Although the modeling of the source extension as a Gaussian in the likelihood is not ideal, it 
is straightforward and computationally fast. The exact shape of the sources is not important 
for small signals; we may be able to analyze the shape with more detail depending on the 
intensity of any signal. 



We perfor med a stacking searc h for five nearby CGs mentioned above following the model 
predictions of iMurase et al.l ( 120081 ) as weights in the likelihood. The size of the clusters in 
the likelihood fit was allowed to vary discretely between the optimal widths for each CR 
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distribution model. The optimal width and differential ffux for each source and each 
model are given in Table [2l The differential fluxes are parametrized as broken power laws: 

^ fTeV-cm-s-^1 - I ^ ' (^/T^V)"- : E < E,,eak .g. 
The parameter B = A ■ E^'^~J} after enforcing continuity at the break energy. 



7. Systematic Errors 



The main systematic uncertainties on the flux limits come from photon propagation 
in ice, absolute DOM sensitivity, and uncertainties in the Earth density profile as well as 
muon energy loss. All numbers are for an neutrino spectrum of muon neutrinos. We 
evaluate the systematic uncertainty due to photon propagation by performing dedicated 
simulat ions with scattering ari d absorption coefficients varied within their uncertainties of 
±10% ( lAckermann et al.l 120061 ). The maximum difference was between the case where both 
scattering and absorption were increased by 10% and the case where both were decreased 
by 10%. The deviation in the observed number of events between these two cases was 11%. 
The range of uncertainty in the DO M sensitivity is take n as ±8%, based on the measured 
uncertainty in the PMT sensitivity (lAbbasi et al.ll2010a( ). Similarly, another dedicated sim- 
ulation where we varied the DOM sensitivity inside the uncertainty leads to a maximum 
systematic uncertainty on the number of events of 9%. These uncertainties on the flux var- 
ied by only about 2% between the northern and southern sky, so only averages over the 
whole sky are reported. Finally, uncertainties in muon energy losses, the neutrino-nucleon 
cross-section, and the rock den sity near the detector in troduce an 8% systematic uncertainty 
for vertically up-going events (lAchterberg et al.l 120071 ). These events are the most affected, 
and this uncertainty is applied to all zeniths to be conservative. A sum in quadrature of 
the systematic uncertainties on the flux gives a total of 16% systematic uncertainty in the 
signal simulation. These systematic uncert ainties are incorpora ted into the upper limit and 
sensit ivity calculations using the method of IConrad et al.l (120031 ) with a modification by iHill 

feoosh . 



The analyses described in Sec. |6] give reliable statistical results (p- values) due to the 
ability to generate background-only data sets by scrambhng the data in right ascension. By 
using the data to estimate background, the systematic errors come only from signal and 
detector simulation used to calculate flux upper limits, sensitivities, and discovery poten- 
tials. Nevertheless, the comparison of background simulation (both atmospheric neutrinos 
and atmospheric muons) with data is an important check to build confidence in the signal 
simulation. 
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Table 2. Galaxy cluster parameters. 



Source 


r.a. [°] 


dec. [°] 


Model 


'^s [°] 


A [TeV-i cm-2 s"!] 


71 


72 


Sbrcak [TeV] 








Model A 


2.0 


1.42 X 10'^'^ 


-2.14 


-4.03 


2.16 X 10^ 


Virgo 


186.63 


12.72 


Model B 
Isobaric 


4.0 
3.0 


1.18 X 10-12 
7.57 X 10"" 


-2.14 
-2.14 


-4.03 
-4.03 


2.16 X 10« 
2.16 X 10'^ 








Central AGN 


0.0 


6.47 X 10-12 


-2.42 


-4.24 


2.13 X 10^ 








Model A 


0.25 


2.78 X 10-" 


-2.14 


-4.03 


2.15 X 10'^ 


Centaurus 


192.20 


-41.31 


Model B 
Isobaric 


0.5 
0.25 


2.20 X 10-" 
1.09 X 10-" 


-2.14 
-2.15 


-4.03 
-4.07 


2.15 X 10'^ 
2.33 X 10^ 








Central AGN 


0.0 


5.10 X 10-" 


-2.45 


-4.28 


2.39 X 10" 








Model A 


0.0 


5.83 X 10-1* 


-2.15 


-4.07 


2.32 X lO'* 


Perseus 


49.95 


41.52 


Model B 
Isobaric 


0.5 
0.0 


4.60 X 10-1* 
6.17 X 10-" 


-2.15 
-2.15 


-4.07 
-4.07 


2.32 X 10" 
2.32 X 10" 








Central AGN 


0.0 


5.97 X 10-" 


-2.40 


-4.20 


1.88 X 10" 








Model A 


0.25 


2.14 X 10-1* 


-2.14 


-4.03 


2.12 X 10" 


Coma 


194.95 


27.94 


Model B 
Isobaric 


0.25 
0.25 


1.34 X 10-1* 
1.83 X 10-" 


-2.14 
-2.15 


-4.03 
-4.07 


2.12 X 10" 
2.30 X 10" 








Central AGN 


0.0 


2.13 X 10-" 


-2.41 


-4.20 


1.89 X 10" 








Model A 


0.0 


4.87 X 10-1* 


-2.15 


-4.07 


2.29 X 10" 


Ophiuchus 


258.11 


-23.36 


Model B 
Isobaric 


0.5 
0.0 


1.50 X 10-1* 
5.50 X 10-" 


-2.15 
-2.15 


-4.07 
-4.11 


2.29 X 10" 
2.49 X lO" 








Central AGN 


0.0 


2.55 X 10-" 


-2.43 


-4.24 


2.12 X 10" 



Note. — as is the optimized sig ma of a 2-dimensional Gaussian distribu tion used in the likelihood. Numerical cal- 
culati ons of the differential fluxes llMurasellpersonal communication, 20091) for each model described in lMurase et al.l 
ll2008h are fit well to broken power laws, parametrized in Eq. |9] 
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Muons created in CR showers in the atmosphere at zenith angles around 80° (near the 
horizon) must travel about 15 km to reach the bottom of IceCube. Only very high energy 
muons can travel such distances. For the simulation to produce the correct zenith distribution 
for these nearly horizontal events, CR composition can be important since protons can 
produce higher energy muons than iron nuclei with the same energy. Two CR models were 
simulated: First, a shghtly rnodifie d version of the poly-gonato model of the galactic CR flux 
and composition (IHoerandell l2003l ) . Above the galactic model cutoff at Z x 4 x 10^^ eV, a 
flux of pure iron is used with an spectrum. As discussed in 13. 2[ these corrections are 
required to reproduce the measured CR spectrum at these energies. Second, a simpler pure 
proton and iron two -component model with a muc h higher contribution of protons is also 
used for comparison ( IGlasstetter fc Hoerandellll999[ ). The final zenith distribution with each 
of these models is shown in Fig. The atmospheric muon simulation is not only affected 
by the primary composition uncertainties at high energy; it is also affected by the hadronic 
model, affecting the production rate of muons at the level of 15% in the region of interest 



for Ic eCube, greater than about 1 TeV, as discussed in the SIBYLL model paper (lAhn et al. 



2009r) and in the compari son between different hadronic models used in CORSIKA presented 



m 



Berghaus et al.l ( l2008l ). For the up-going region, several models of atmospheric neutrino 



fluxes, both conventional fluxes from pion and kaon decay and prompt fluxes from charmed 
meson decay, are shown in Fig. [161 To represe nt the low and high predictions, conventional 
and prompt models a re used in pairs: H onda ( iHonda et al.l 120071 ) for the conventional flux 
paired with Sarcevic (Enberg et al.ll2008[ ) for the prompt flux represent the l ow prediction , 
and Bartol (IBarr et al.l 120041 ) for the conventional flux paired with Naumov ( lBugaevlll989[ ) 
for the prompt flux represent the high prediction. Additional uncer tainty in the pre dicted 
atmospheric neutrino rate is estimated to be about 40% at 1 TeV ( iBarr et al.ll2006l ). We 
conclude that our data agree with background simulation at the final level, within the range 
of uncertainties allowable by existing CR composition and atmospheric neutrino models. 



8. Results 

The results of the all-sky scan are shown in the map of the pre-trial p- values in Fig. [T71 
The most significant deviation from background is located at 113.75° r.a., 15.15° dec. The 
best-fit parameters are Ug = 11.0 signal events above background, with spectral index 7 = 
2.1. The pre-trial estimated p- value of the maximum log likelihood ratio at this location is 
5.2 X 10^^. In trials using data sets scrambled in right ascension, 1,817 out of 10,000 have 
an equal or higher significance somewhere in the sky, resulting in the post-trial p-value of 
18%. Upper limits for an E'"^ flux of z/^ + are presented in Fig. [181 In all cases, an equal 
flux of neutrinos and anti-neutrinos is assumed. 



- 32 - 



w4000 
I35OO 

LU 

3000 
2500 
2000 
1500 
1000 
500 




Data 

Atm fi MC (poly-gonato + iron) 
Atm fi MC (2-component) 
^^^H Atm V MC (Bartol+Naumov) 
Atm V MC (Honda+Sarcevic) 
I I Aim MC (low models) 
I I Atm n+v MC (high models) 




cos(e) 



Fig. 16. — Distribution of reconstructed cosine zenith for the final event sample compared to 
the models discussed in the text. Honda and Sarcevic are summed with poly-gonato to rep- 
resent the set of low predictions, and Bartol and Naumov are summed with the 2-component 
model for the high predictions. Only statistical errors are shown. The two-component model 
has limited statistics, causing the peaks and valleys. Systematic un certainties of neu trino 
production in CR showers are estimated to be a bout 40% at 1 TeV fjBarr et al.ll2006l) and 
15% in the muon rate greater than about 1 TeV (lAhn et al.ll2009l : iBerghaus et al.ll2008l ). 



The results of the point-source search in the direction of 39 source candidates selected 
a priori are given in Tableland also shown in Fig. [19] with the IceCube median sensitivity. 
The most significant source on the list was PKS 1622-297 with a pre-trial estimated p-value 
of 5%. The post-trial p-value of 62% was again determined as the fraction of scrambled 
data sets with at least one source with an equal or higher significance. The mean number 
of events at the final cut level required for the discovery of a point source is also shown in 
Fig. EHl along with the average background in a circular bin wit h 1° radius. Included in 



Fig. [I9]is a preliminary comparison to the ANTARES experiment ( ICoylell2010l ). ANTARES 
is primarily sensitive to GeV-TeV energy neutrinos in the southern sky, so the coverage in 
energy is quite complementary to this IceCube analysis. 
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Fig. 17. — Equatorial skymap (J2000) of pre-trial significances (p- value) of the all-sky point 
source scan. The galactic plane is shown as the solid black curve. 




Fig. 18. — Equatorial skymap (J2000) of upper limits of Feldman- Cousins 90% confidence 
intervals for an E~'^ flux of z/^ + z/^. The galactic plane is shown as the solid black curve. 
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Table 3. Results for the a priori source candidate list. 



01).ieet 


r.a. n 


dec. [°] 


^90 




])- value 






Nr. 


B,- 


Cyg 0B2 


308.08 


41.51 


6.04 


10.54 




0.0 




2 


1.8 


MGRO J2019+37 


305.22 


36.83 


7.50 


13.3 


0.44 


1.0 


2.8 


2 


1.9 


MGRO J1908+06 


286.98 


6.27 


3.73 


6.82 


0.43 


1.5 


3.9 


4 


3.1 


Gas A 


350.85 


58.81 


9.04 


15.92 




0.0 




1 


1.8 


IG443 


94.18 


22.53 


3.80 


6.62 




0.0 




1 


2.0 


Gcminga 


98.48 


17.77 


3.91 


6.66 


0.48 


0.7 


2.1 


1 


2.3 


1 ATI 1 

Crab iNcbma 


83.63 


22.01 


3.70 


6.58 




0.0 




1 


2.0 


lES 1959+650 


300.00 


65.15 


10.74 


19.18 




0.0 







2.0 


lES 2344+514 


356.77 


51.70 


7.24 


12.96 




0.0 







1.8 


3C66A 


35.67 


43.04 


10.89 


19.70 


0.24 


3.4 


3.9 


3 


1.9 


H 1426+428 


217.14 


42.67 


6.14 


10.94 




0.0 




3 


1.8 


BL Lac 


330.68 


42.28 


10.80 


18.70 


0.25 


2.6 


3.9 


3 


1.8 


Mrk 501 


253.47 


39.76 


8.11 


14.14 


0.41 


1.3 


3.9 


3 


2.0 


TV T 1 A C\ 

Mrk 421 


166.11 


38.21 


11.71 


20.14 


0.15 


2.6 


1.9 


2 


2.0 


W Comae 


185.38 


28.23 


4.46 


8.06 




0.0 







1.9 


lES 0229+200 


38.20 


20.29 


6.89 


12.06 


0.19 


4.0 


2.8 


4 


2.1 


M87 


187.71 


12.39 


3.42 


5.98 




0.0 




2 


2.5 


S5 0716+71 


110.47 


71.34 


13.28 


23.56 




0.0 







1.6 


M82 


148.97 


69.68 


19.14 


32.84 


0.4 


2.0 


3.9 


4 


1.8 


3G 123.0 


69.27 


29.67 


5.59 


10.66 


0.44 


1.3 


2.7 


1 


1.9 


3G 454.3 


343.49 


16.15 


3.42 


5.92 




0.0 




1 


2.3 


4G 38.41 


248.81 


38.13 


6.77 


11.86 


0.48 


0.9 


3.9 


2 


2.0 


PKS 0235+164 


39.66 


16.62 


6.77 


11.62 


0.15 


5.3 


3.0 


5 


2.3 


PKS 0528+134 


82.73 


13.53 


3.63 


6.72 




0.0 




2 


2.4 


PKS 1502+106 


226.10 


10.49 


3.26 


5.78 




0.0 







2.5 


3C 273 


187.28 


2.05 


3.61 


6.54 




0.0 




3 


3.4 


NGC 1275 


49.95 


41.51 


6.04 


10.54 




0.0 




2 


1.8 


Cyg A 


299.87 


40.73 


7.84 


13.44 


0.46 


1.0 


3.5 


3 


1.9 


IG-22 maximum 


153.38 


11.38 


3.26 


5.86 




0.0 




1 


2.5 


Sgr A* 


266.42 


-29.01 


80.56 


139.26 


0.41 


1.1 


2.7 


4 


3.3 


PKS 0537-441 


84.71 


-44.09 


113.90 


201.82 




0.0 




3 


3.5 
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22 Strings Sensitivity 275.7 d 

40 Strings Discovery Potential 375.5 d 




sin (8) 



Fig. 19. — Median sensitivity to a point-source flux of t'^+z/^ as a function of declinati on, 
shown for the 22-string IceCube southern and northern sky analyses (lAbbasi et al.ll2009al lbl). 
this 40-string analysis, and preliminary estimated sensitivities fo r 1 year for t he ANTARES 
experiment, primarily sensitive in the GeV-TeV energy range, ( ICoyld |2010| ) and the final 
IceCube configuration (using the event selection based on this work for the up-going region). 
For the a priori source list, upper limits of Feldman-Cousins 90% confidence intervals for an 
E~'^ flux of + Ufj, are shown. In addition, we show the discovery potential for this work. 



The Milagro TeV source stacking search resulted in a post-trial p-value of 32% with 
best fit rig = 7.6 and spectral index 7 = 2.6. The starburst galaxy stacking search resulted 
in an under-fiuctuation with best fit = 0, since we do not allow negative values of Ug. 
Finally, the CGs stacking search yielded a post-trial p-value of 78% with Us = 1.8. These 
results, along with sensitivities and upper limits, are summarized in Table HI 



9. Implications for Models of Astrophysical Neutrino Emission 



The IceCube Neutrino Observatory aims to further our understanding of astrophysical 

phenomena, constraining models even in the absenc e of a detectio n. Fig. [2T] shows our 

sensitivity to some specific predictions. The model of iMorlino et al.l (120091 ) is for SNR RX 
J1713. 7-3946. This analysis is largely insensitive to spectra which cut off below 100 TeV 
in the southern sky. Applying this emission model at the location of the Crab Nebula 
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Table 3 — Continued 



Object 


r.a. [°] 


dec. [°] 


$90 


^90 


p- value 


ns 


7 


Nio 


Bio 


Cen A 


201.37 


-43.02 


109.51 


191.56 




0.0 




4 


3.5 


PKS 1454-354 


224.36 


-35.65 


92.56 


156.74 




0.0 




4 


3.5 


PKS 2155-304 


329.72 


-30.23 


105.41 


182.90 


0.28 


1.7 


3.9 


3 


3.4 


PKS 1622-297 


246.53 


-29.86 


152.28 


263.86 


0.048 


3.0 


2.6 


4 


3.3 


QSO 1730-130 


263.26 


-13.08 


24.83 


43.30 




0.0 




4 


3.5 


PKS 1406-076 


212.24 


-7.87 


16.04 


28.72 


0.42 


1.3 


2.3 


4 


3.3 


QSO 2022-077 


306.42 


-7.64 


12.18 


21.78 




0.0 




2 


3.3 


3C279 


194.05 


-5.79 


11.94 


21.36 


0.33 


3.6 


3.0 


7 


3.5 



Note. — and are the upper limits of the Feldman- Cousins 90% confi- 

dence intervals for an E~'^ flux normalization of z/^ and + I'r, respectively, in units of 
10-^2 TeV"^cm-2s-i (i.e. dN/dE < ■ (E/TeY)-^). is the best-fit number of signal 
events; when Us > the (pre-trial) p-value is also calculated and the spectral index 7 is 
given. Nio is the actual number of events observed in a bin of radius 1°. The background 
event density at the source declination is indicated by the mean number of background 
events Bio expected in a bin of radius 1°. 



Table 4. Results for the stacked source searches. 



Catalog 


N Sources 


Model 


p- value 


Sensitivity 


Upper Limit 


i^fj, + i^T Sensitivity 


1^1^ + i^T Upper Limit 


Milagro Sources 


17 
6 


E^^, Uniform 
6 SNR Assoc.'' 


0.32 

C 


$90 = 9.0 


4.90 = 12.3 


#9" = 15.8 
SF = 2.9 


*90 = 24.5 
SF = 7.2 


Starburst Galaxies 


127 


E-'-^,oc FIR Flux 




= 33.1 


$9" = 33.1 


#90 = 58.6 


$90 ^ 5g g 






Model A'" 








SF = 8.4 


SF = 7.8 


Clusters of Galaxies 


5 


Model B*^ 
Isobaric'' 
Central AGN^ 


0.78 






SF = 14.4 
SF = 13.2 
SF = 6.0 


SF = 12.0 
SF = 13.2 
SF = 6.0 



Note. — Median sensitivities and upper limits at 90% CL for i/^ and + fluxes are given in two ways: as $90 fQj. £'—2 spectrum, i.e. the total flux 
from all sources dN/dE < <I>90 IQ—^^ TeV~^cm~'^s~^(i?/TeV)~^, or as a scaling factor (SF) relative to the models given in the footnotes. For example, 
if the Central AGN model flux normalization were 6.0 times higher, we would rule it out at 90% CL. All models predict equal fluxes of tau and muon 
neutrinos. 

HHalzen et all l l2008h 

iMurase et"al] ||2008|) . see Table [2] 

'^We did not calculate an a priori p-value for just the six SNR associations discussed in lHalzen et al. I ll2008h . since they are included in the search over all 
Milagro sources. However, some of these sources were the most significant on the list. Analyzed as a single subgroup, an a posteriori p-value of 0.02 was 
found with best fit parameters ris = 15.2 and 7 = 2.9. The true trial factor is incalculable since this was done after unblinding, but these remain sources 
of interest for future searches. 
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Fig. 20. — Mean number of z/^ + z/^ signal events at the final cut level required for a 
discovery at 5cr in 50% of trials and the mean number of background events in a circular bin 
with a radius of 1° versus sine of the declination. 



{6 = 22.01°) we obtain an upper limit that rules out a flux 3.2 times higher than the 
prediction. 

The Mil agro hot spot MGRO J1852+01 is the brightest source of six SNR associations 
considered in lHalzen et al.l (2008). The stacking results were already shown in Table HI Our 
upper limit for just this one brightest source is a factor of 7.9 away from excluding this model 
at 90% CL. The best fit for MGRO J1852+01 is to 7.0 events with 7 = 2.9, which increases 
the upper hmit compared to the average background-only case. 

The nearest AGN, Centaurus A (Gen A), has been discussed as a potential source of 
ultra high energy GRs in the context of results from the Pierre Auger Observatory (PAO). 
The point source likelihood fit at the location of Gen A resulted in zero signal events in 
this analysis, setting an upper limit that is 5.3 times higher than the z/^ prediction by 
Koers fc Tinyakovl (120081 ) for the most optimistic case where the protons have a spectral 



index = 3. 



Star burst galaxies were already presented as sources of interest in Sec. 16.41 Recent de- 
tections (lAbdo et al.l I2OIOI : lAcciari et al.l l2009l : lAharonian et al.l |2005|) of very high energy 
photons from the nearest luminous starburst galaxies M82 and NGG253, each characterized 
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Fig. 21. — Differential ffux of three theoretical models shown with the IceCube 40-string 
up per limit (90% CL) a nd discovery potential in each case. Shown are the i/^ + predictions 
of iMorlino et all (120091) for SNR RX J1713.7-3946 but moved to the location o f the Crab 



Nebula, Halzen et al. 



(120081 ) for MGRO J1852+01, and lKoers fc Tinvakovl feOOSf ) for Cen A 



under the most optimistic condition, where the protons have a spectral index = 3. 



by star-forming regions with high supernova rates in the core, support the belief that the 
observed enhanced gamma ray emission is due to CR interactions. Under the assumption 
that the GeV-TeV photons originate from the decay of neutral pions produced when pro- 
tons that are shock-accelerated by SNRs in the starburst core inelastically scatter against 



target hydrogen atom s with densities of the order of 10^ cm ^ (Ide Cea del Pozo et al.l 12009 



Persic. M. et al.ll2008l). an order-of- magnitude calculation of the resulting flux of muon neu- 
trinos based on lKelner et al.l (120061 ) can be made. The muon neutrino upper limit from M82 
is about 400 times higher than the rough prediction. For NGC 253 in the southern sky, the 
muon neutrino upper limit is about 6000 times higher than the prediction. 



10. Conclusions 

A search for sources of high energy neutrinos has been performed using data taken 
during 2008-09 with the 40-string configuration of the IceCube Neutrino Observatory. Five 
searches were performed: 1) a scan of the entire sky for point sources, 2) a predefined list 
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of 39 interesting source candidates, 3) stacking 16 sources of TeV gamma rays observed by 
Milagro and Fermi, along with an unconfirmed hot spot (17 total sources), 4) stacking 127 
starburst galaxies, and 5) stacking five nearby CGs, testing four different models for the 
CR distribution. The most significant result of the five searches came from the all-sky scan 
with a p-value of 18%. The cumulative binomial probability of obtaining at least one result 
of this significance or higher in five searches is 63%. This result is consistent with the null 
hypothesis of background only. 

The sensitivity of this search using 375.5 days of 40-string data already improves upon 
previous point-source searches in the TcV-PeV energy range by at least a factor of two, 
depending on declination. The searches were performed using a data set of up-going at- 
mospheric neutrinos (northern sky) and higher energy down-going muons (southern sky) 
in a unified manner. During 2010-2011, 79 strings of IceCube are operating and detector 
construction should finish during the austral summer of the same years. The full IceCube 
detector should improve existing limits by at least another factor of two with one year of 
operation. Additional improvement is foreseeable in the down-going region by developing 
sophisticated veto techniques and at lower energies by using the new dense sub-array, Deep- 
Core, to its fullest potential. 
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